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We present a comprehensive study of event-by-event multiplicity fluctuations in nucleon-nucleon 
and nucleus-nucleus interactions from AGS/FAIR to RHIC energies within the UrQMD transport 
approach. The scaled variances of negative, positive, and all charged hadrons are analysed. The 
scaled variance in central Pb-|-Pb collisions increases with energy and behaves similar to inelastic 
p-|-p interactions. We find a non-trivial dependence of multiplicity fluctuations on the rapidity 
and transverse momentum interval used for the analysis and on the centrality selection procedure. 
Quantitative predictions for the NA49 experiment are given, taking into account the acceptance of 
the detector and the selection procedure of central events. 
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I. INTRODUCTION 

At high energy densities 1 GeV/fm'^) a phase tran- 
sition from a hadron gas to a quark-gluon-plasma (QGP) 
is expected to occur. There are indications that at RHIC 
and top SPS energies a quark-gluon-plasma is created in 
the early stages of heavy ion collisions P, And in- 
deed, the energy dependence of various observables show 
anomalies at low SPS energies which might be related to 
the onset of deconfinement [1, [3. 

While several observables Q have been proposed 
throughout the last decades to study the characteristics 
of the highly excited matter created in heavy ions col- 
lisions the ones related to fluctuations and correlations 
seem to be the most prospective. Fluctuation probes 
might be more adequate for the exploration of heavy 
ion reactions, because the distributions of energy den- 
sity or initial temperature, isospin and particle density 
have strong fluctuations from event to event 0, 0. 
On the theoretical side event-by-event fluctuations were 
suggested to study 

• kinetic and chemical equilibration in nuclear colli- 
sions [i, M, M, E [H, El, m, [M [13 > 

• the onset of the deconfinement phase [1, [H, [l^ [20l . 

• the location of the tri-critical end-point of the QCD 
phase transition [H, [13, [111 or 

• the formation of exotic states, like DCCs 12911. 
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On the experimental side, progress has been made by 
many experiments to extract momentum and particle 
number ratio fluctuations from heavy ion reaction: Cur- 
rently, event-by-event fluctuations are actively studied in 
the SPS energy regime (star ting from 20A GeV on) by 
the NA49_group [30^ [H, [M [M [H, [H, [H [H , the 
CERES [3a[40,[4l[l4fand the WA98 collaboration [H. 
At RHIC energies the PHENIX [H, [H, HI and STAR 
(47I list |49| experiments are addressing the field of single 
event physics. 

In [2^ it was predicted that the onset of deconfine- 
ment should lead to a non-monotonous behaviour in mul- 
tiplicity fluctuations ("shark fin"). Also droplet forma- 
tion during the phase transition is expected to produce 
non-statistical fluctuations 10-100 times the Poisson ex- 
pectation [50l|. Furthermore, lattice QCD calculations 
suggest the existence of a critical point in the phase dia- 
gram of strongly interacting matter which separates the 
first order phase transition from a crossover. Thus, if the 
system passes the vicinity of the critical region during 
its evolution and remains there for long enough time one 
expects an increase of multiplicity fluctuations ^27,] . 

The NA49 collaboration is currently searching for such 
anomalies in the energy dependence of multiplicity fluc- 
tuations in Pb-|-Pb collisions. A similar program to 
search for the critical point and signals for the onset of 
deconfinement will be undertaken by the RHIC exper- 
iments (the planned critRHIC program ) with a lower- 
ing of the RHIC's beam energy towards the SPS energy 
regime and the NA61 (SHINE) experiment [5l| with the 
focus on light ion collisions. For the present investigation 
however, we will focus on the soon available data from 
the NA49 experiment on multiplicity fluctuations. Un- 
fortunately both the geometrical acceptance of the detec- 
tor and the centrality selection in the NA49 experiment 
is not trivial and have an influence on the multiplicity 
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fluctuations. In order to observe an increase of fluctu- 
ations caused by one of the effects mentioned above, a 
systematic theoretical investigation within a transport 
approach is needed. Only with this baseline for the ex- 
pected multiplicity fluctuations within the experimental 
acceptance a possible excess of fluctuations in data could 
be unambiguously interpreted as a signal for the critical 
point or the onset of deconfinement. The model predic- 
tions presented in this paper are obtained using UrQMD 
1.3 [i^,!!^- For a complementary transport theoretical 
study of multiplicity fluctuations, the reader is referred 
to [5|,[5i,[5i,[53. 

II. MEASURE OF MULTIPLICITY 
FLUCTUATIONS 

The probability to have in an event a given number 
of particles n in the acceptance is denoted as P{n), with 
the normalisation "^Pin) — 1. 

n 

The measure of multiplicity fluctuations used in this 
paper is the scaled variance lu defined as 

Varju) (n^) - {n)' 
(n) (n) 

where Var{n) — — {n))'^P(n) and (n) —^n- P{n) 

n n 

are the variance and the mean of the multiplicity distri- 
butions, respectively. 

This measure is used because of its two properties. 
Firstly in a grand-canonical statistical model neglecting 
quantum effects the multiplicity is a Poisson distribution: 

P(n) = M!.exp(-(n)) (2) 

The variance of a Poisson distribution is equal to its 
mean, the scaled variance is therefore to — 1, indepen- 
dent of mean multiplicity. 

Secondly in a Wounded Nucleon Model [S^ , the scaled 
variance in A + A collisions is the same as in proton- 
proton collisions provided the number of wounded nucle- 
ons is fixed. If the particles are produced independently 
in momentum space, the scaled variance in a limited ac- 
ceptance is related to the scaled variance in full phase 
space (47r): 

l^acc = P ■ (l^47r - 1) + 1 (3) 

where p is the fraction of tracks which are in the corre- 
sponding acceptance. For a small acceptance p the scaled 
variance approaches 1. Note that effects like resonance 
decays, quantum statistics and energy-momentum con- 
servation introduce correlations in momentum space and 
therefore a scaling according to equation [3] is generally 
not valid. 

In the following the scaled variance of positively, neg- 
atively and all charged hadrons are denoted as uj{h^), 
uj{h~) and u!{h^), respectively. 



III. THE URQMD MODEL 

For our investigation, we apply the Ultra-relativistic 
Quantum Molecular Dynamics model (UrQMD vl.3) [H, 
to heavy ion reactions from i^boam = 20A GeV to 
£^bcam = 158A GeV. This microscopic transport ap- 
proach is based on the covariant propagation of con- 
stituent quarks and di-quarks accompanied by mesonic 
and baryonic degrees of freedom. It simulates multiple 
interactions of in-going and newly produced particles, the 
excitation and fragmentation of colour strings and the 
formation and decay of hadronic resonances. Towards 
higher energies, the treatment of sub-hadronic degrees of 
freedom is of major importance. In the present model, 
these degrees of freedom enter via the introduction of a 
formation time for hadrons produced in the fragmenta- 
tion of strings [sl, [13, [6l| . A phase transition to a quark- 
gluon state is not incorporated explicitly into the model 
dynamics. However, a detailed analysis of the model in 
equilibrium, yields an effective equation of state of Hagc- 
dorn type jH, [gl] . 

This model has been used before to studye vent-by- 
event fluctuations rather successfully [1, [ill, US [IE [H, 
[stI . [6^ . [65| and yields a reasonable description of inclu- 
sive particle distributions. For a complete review of the 
model, the reader is referred to (s^. [53|. 



IV. ENERGY DEPENDENCE OF 
MULTIPLICITY FLUCTUATIONS 

The energy dependence of the mean multiplicity, nor- 
malised by the number of nucleons of one projectile 
{A = 1 for p-fp, p+n, A = 208 for Pb-hPb) of posi- 
tively, negatively and all charged particles in p-l-p, p-l-n 
and Pb-I-Pb collisions is shown in flgure[T] For p-l-p and 
p-|-n interactions all inelastic collisions are selected. For 
Pb+Pb the impact parameter of the collisions are set 
to 6 = 0. The calculations were performed for AGS 
{Eiab = 6.87A GeV), SPS {Eiab = 20^, 30/1, 40^, 80^ 
and 158^ GeV) and RHIC (^/i^ = 62.4 and 200 GeV) 
energies. In the UrQMD 1.3 model the mean multiplicity 
per number of projectile nucleons is significantly larger 
in Pb-|-Pb collisions in comparison to p-|-n interactions. 

The mean multiplicity of all charged hadrons in p-fp 
interactions obtained by various experiments is parame- 
terised in [6^ as 

< n± >w -4.2 + 4.69 • {^/^/GeVf-^^ (4) 

Except for top RHIC energies the parametrisation of the 
experimental data is in agreement with the UrQMD re- 
sult. 

The energy dependence of scaled variance in full phase 
space is shown in figure [2] An increase of scaled vari- 
ance with increasing collision energy is observed for p-l-p, 
p-l-n and Pb-|-Pb collisions. For AGS and low SPS en- 
ergies the scaled variance is smaller than one and the 
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FIG. 1: (Color online) Mean multiplicity in Att of inelastic 
p+p, p+n and central Pb+Pb collisions as a function of col- 
lision energy. Top left: positively, top right: negatively, bot- 
tom: all charged hadrons. 

multiplicity distributions are narrower than the corre- 
sponding Poisson distributions. For higher energies the 
scaled variance is larger than one. A similar behaviour 
of p+p and p-fn collisions is observed, the small differ- 
ence is probably caused by the additional proton in p-l-p 
collisions, which does not fluctuate. Therefore the scaled 
variance for positively and all charged particles is a bit 
lower in p-l-p than in p-t-n collisions. The scaled variance 
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^ (GeV) 




^ (GeV) 

FIG. 2: (Color online) Scaled variance in Att of inelastic p-l-p, 
p+n and central Pb+Pb collisions as a function of collision 
energy in comparison to hadron gas model predictions [g^] 
for Pb+Pb collisions. Top: positively, middle: negatively, 
bottom: all charged hadrons. 



in Pb+Pb collisions behaves similar as in p+P interac- 
tions. The HSD model yields similar results [55| . 

For positively and negatively charged hadrons the 
scaled variance is similar, where the values are about 
twice as high for all charged hadrons. This is partly due 
to resonances decaying into two oppositely charged par- 
ticles. Such a resonance is detected as two charged par- 
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tides, therefore the fluctuations are increased. 

The experimental data on the energy dependence of 
scaled variance of all charged hadrons in p+p interactions 
is parametrised in [66| as 



« 0.35 • 



(< > -1)^ 
< n± > 



(5) 



The UrQMD results are in agreement with the data at 
AGS and SPS energies, but the fluctuations are slightly 
overpredicted at RHIC energies. 

Figure [2] shows that the energy dependence of scaled 
variance predicted by the UrQMD model is totally dif- 
ferent to the predictions of a hadron gas model 16711 (for 
further details the reader is referred also to [68l. l69l. ItoI . 
[ZH, III)- In the grand-canonical (GCE), canonical (CE) 
and micro-canonical (MCE) ensemble the scaled variance 
stays constant for high energies, where in the UrQMD 
model it strongly increases with energy. Therefore ex- 
perimental data on multiplicity fluctuations, preferably 
at high (RHIC, LHC) energies, should be able to dis- 
tinguish between hadron gas and string-hadronic mod- 
els |55| . 

For a more differential study of fluctuations and for a 
better comparison to experimental results, three differ- 
ent rapidity intervals, one at midrapidity < ?; < 1, one 
at forward rapidity 1 < y < ybeam and a combination 
of both < y < ybeam, covering most of the forward 
hemisphere, were taken. For p + n collisions the for- 
ward hemisphere includes the rapidity of the projectile 
neutron. The scaled variance for these intervals for pos- 
itively, negatively and all charged particles are shown in 
figures [3][5l 

As in full phase space, in the three different rapidity 
intervals a similar behaviour of scaled variance of p -|- p, 
p + n and Pb-f Pb collisions was observed. 

The energy dependence of ffuctuations in the forward 
hemisphere (0 < y < ybeam) looks similar to the one in 
the full phase space, the absolute number of scaled vari- 
ance is similar to the result expected when applying the 
acceptance extrapolation according to formula [3] (shown 
as stars in figures . 

For low energies a large fraction of particles is in the 
midrapidity interval {0 < y < 1) where a very small 
amount of particles is in the forward rapidity interval 
{1 < y < ybeam)- With increasing energy both the width 
and the number of particles in the forward rapidity in- 
terval increases strongly, where the number of particles 
in the forward rapidity interval increases only weakly. 

At midrapidity (0 < y < 1) the scaled variance is in 
the same order of magnitude as in the rapidity inter- 
val < y < ybeam, but the mean multiplicity is much 
lower. The acceptance extrapolation formula [3] strongly 
underpredicts fluctuations in this rapidity region. At for- 
ward rapidity (1 < y < ybeam) the fluctuations are much 
smaller than predicted by the acceptance extrapolation 
formula. For lower energies the scaled variance decreases 
with energy, for higher energies it increases. This can be 
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FIG. 3: (Color online) Scaled variance of positively charged 
hadrons produced in inelastic p + p, p + n and central 
Pb+Pb collisions as a function of collision energy. Top: 
< y < ybeam, middle: < y < 1, bottom: 1 < y < ybeam- 



qualitatively understood by the interplay of an increas- 
ing fraction of particles in this rapidity interval and an 
increasing scaled variance in 47r, which is smaller than 1 
for lower and larger than 1 for higher energies. 



5 



3 
2.5 

2 
1.5 

1- 

0.5 


1 

3 
2.5 

2 
1.5 

1 

0.5 



-•— Pb+Pb, b=0 

★ Pb+Pb 4jt, acc. extrap. 

▲ P+P 

▼ p+n 




3 
2.5 

2 
1.5 

1 

0.5 



1 



10 



10 



10^ 
^ (GeV) 



— •— Pb+Pb, b=0 
★ Pb+Pb 4ii, acc. extrap 
▲ P+P 
▼ p+n 


* 







10 


10=^ 


(GeV) 


— •— Pb+Pb, b=0 
★ Pb+Pb 4it, acc. extrap. 
▲ P+P 
▼ p+n 




★ 








^ — 





10^ 
^ (GeV) 



FIG. 4; (Color online) Scaled variance of negatively charged 
hadrons produced in inelastic p + p, p + n and central 
Pb+Pb collisions as a function of collision energy. Top; 
< y < ybeam, middle: < y < 1, bottom: 1 < j/ < i/beam- 
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FIG. 5: (Color online) Scaled variance of all charged hadrons 
produced in inelastic p+p, p + n and central Pb+Pb collisions 
as a function of collision energy. Top: < y < ybeam, middle: 
< y < 1, bottom: 1 < y < ybeam- 



V. RAPIDITY AND TRANSVERSE 
MOMENTUM DEPENDENCE 

As already showed in section IIVI the scaled variance 
is a non-trivial function of the selected phase-space. In 
order to study the dependence of scaled variance on ra- 
pidity, 12 different rapidity intervals are constructed in 



such a way that the mean multiplicity in each interval is 
the same. If the scaled variance would follow the accep- 
tance scaling formula[ni the scaled variance would be the 
same in each interval. In figure [6] it is shown that this 
is not the case. The scaled variance is much higher near 
midrapidity than in forward and backward rapidities. 

The transverse momentum dependence of scaled vari- 
ance is shown in figure [7] for the full longitudinal phase- 
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FIG. 6: Rapidity dependence of scaled variance in UrQMD 
simulation performed in full acceptance of positive (top) , neg- 
ative (middle) and all charged (bottom) hadrons in central 
Pb+Pb collisions at 158A GeV . The rapidity bins are con- 
structed in such a way that the mean multiplicity in each bin 
is the same. 



space and for a midrapidity and a forward rapidity inter- 
val. The scaled variance decreases with increasing trans- 
verse momentum for the full acceptance and at forward 
rapidity. At midrapidity it stays approximately constant. 
The decrease of scaled variance is stronger for positively 
charged hadrons than for negatively charged ones be- 



cause the protons, which have smaller relative fluctua- 
tions due to the large number of protons which enter the 
collision, have a larger mean transverse momentum. 

A similar effect of decreasing fluctuations for larger 
rapidities and transverse momenta is observed as a result 
of energy- and momentum conservation in a hadron gas 
model using the micro-canonical ensemble fl^. It costs 
more energy to create a particle with high momentum, 
therefore their number is expected to fluctuate less. 



VI. PREDICTIONS FOR THE NA49 
EXPERIMENT 

Preliminary data of the NA49 experiment on the en- 
ergy dependence of multiplicity fluctuations in very cen- 
tral Pb-|-Pb collisions was shown in [tI, [T^ . Final data 
obtained in a larger geometrical acceptance will be pub- 
lished soon. 

In order to compare the experimental data with model 
calculations, both the geometrical acceptance of the de- 
tector and the centrality selection have to be imple- 
mented in the model calculation. 

The geometrical acceptance of the NA49 experi- 
ment [76| is located mostly in the forward hemi- 
sphere. The acceptance defined by the detec- 
tor geometry and the track selection criteria is dif- 
ferent for each collision energy and a complicated 
function of the particle momentum p. Accep- 
tance tables in y(7r), px and (f> can be obtained 
at the author's website (http://www.ikf.physik.uni- 
frankfurt.de/users/lungwitz/acceptance/). In the NA49 
experiment it is not possible to identify a particle on the 
track-by-track basis, therefore for the calculation of ra- 
pidity in the fixed target laboratory system pion mass is 
assumed. The rapidity is then transformed into the cen- 
ter of mass system of the collision. 

For the UrQMD model predictions showed in this sec- 
tion both the geometrical acceptance defined by the ac- 
ceptance tables and the assumption of pion mass when 
calculating rapidity and transforming into the center of 
mass system are taken into account. 

In the NA49 experiment the centrality of a collision can 
be measured by the energy of projectile spectators, which 
are registered by a calorimeter. This veto calorimeter is 
adjusted in such a way that it registers all spectator pro- 
tons, neutrons and fragments of the projectile. The lower 
the energy in the veto calorimeter the more central is the 
collision. For the multiplicity fiuctuation analysis the 1% 
most central collisions are selected. For this selection, the 
fluctuations in the number of target participants is also 
minimized 57]. Remarks on the contribution of target 
participant fluctuations to multiplicity fluctuations are 
presented in (77| . 

A small fraction of the produced particles in 
a collision is also entering the calorimeter and 
introducing a small bias on centrality measure- 
ment. Acceptance tables of the veto calorimeter 
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FIG. 7: (Color online) Transverse momentum dependence of 
multiplicity fluctuations of positively (top), negatively (mid- 
dle) and all charged hadrons (bottom) for all rapidities), 
< y < 0.5 and 1.25 < y < 1.75. The transverse momentum 
bins are constructed in such a way that the mean multiplicity 
in each bin is the same. 



as a function of p, pT and can be obtained 
on the author's website (http://www.ikf.physik.uni- 
frankfurt.de/users/lungwitz/acceptance/) and are used 
for the UrQMD predictions of multiplicity fluctuations. 

The predictions for the energy dependence of scaled 
variance measured in the NA49 experiment are shown in 



FIG. 8: UrQMD predictions for scaled variance of positively 
charged hadrons in NA49 acceptance produced in central 
Pb-fPb collisions as a function of collision energy. Top: 
< y < Vbeam, middle: < y < 1, bottom: 1 < y < ybeam- 



figures [SlflOl To study the influence of centrality selection 
the scaled variance was also calculated for Pb-|-Pb colli- 
sions with a zero impact parameter b. 

The UrQMD model predicts a weak energy dependence 
of scaled variance for positively and negatively charged 
hadrons in forward acceptance. At midrapidity, at full 
experimental acceptance and for all charged hadrons at 
all acceptances an increase of scaled variance with colli- 
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FIG. 9: UrQMD predictions for scaled variance of nega- 
tively charged hadrons in NA49 acceptance produced in cen- 
tral Pb-|-Pb collisions £is a function of collision energy. Top: 
<y < ybeam, middle: < y < 1, bottom: 1 < y < ybeam- 
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FIG. 10: UrQMD predictions for scaled variance of all charged 
hadrons in NA49 acceptance produced in central Pb+Pb col- 
lisions as a function of collision energy. Top: < y < ybeam, 
middle: < y < 1, bottom: 1 < y < ybeam- 



sion energy is predicted. 

In forward acceptance a UrQMD simulation for events 
with zero impact parameter {h = 0) gives similar results 
to the simiilation for events selected according to their 
energy in the veto calorimeter. In midrapidity and in full 
experimental acceptance the scaled variance for events 
selected by their veto energy is larger, probably due to 
target participant fluctuations. 



VII. SUMMARY 

We present predictions for evcnt-by-evcnt multi- 
plicity fluctuations for very central Pb-|-Pb and nu- 
cleon+nucleon interactions from Eiab = 20A GeV to 
Eiab = 158^ GeV within a hadron-string transport ap- 
proach. We find that the fluctuations do generally in- 
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crease strongly towards higher beam energies, both for 
elementary and massive nuclear reactions. This can be 
used to distinguish the present model from a hadron gas 
model, which predicts a rather weak dependence of the 
fluctuations on energy. The amount of fluctuations in full 
phase space is generally non-Poissonian (smaller at low 
energies, higher at high energies), crossing the Poissonian 
value in the SPS energy regime. Applying a forward ra- 
pidity cut, yields a non-monotonous behaviour as a func- 
tion of energy, with a local minimum at SPS energies. 
We further predict the rapidity dependence of the scaled 
variance at the highest SPS energy and find a strong 
rapidity dependence of the fluctuations, even if trivial 
multiplicity effects are scaled out. This might render the 
procedure to simply scale thermal model predictions for 
fluctuations to the experimentally covered phase space 
questionable. The transverse momentum dependence of 
the fluctuations does generally tend to decrease towards 
higher transverse momenta. This effect is related to en- 
ergy conservation, allowing stronger fluctuations for low 
energetic particles, while constraining the high energetic 
particles. Finally, we analyse the influence of the veto 
trigger used by the NA49 experiment, compared to sim- 
ple zero impact parameter reactions. Here we observe a 
systematic deviation for the midrapidity results on the 
order of 10%, at forward rapidities, the veto trigger can 
be well approximated with a the zero impact parameter 
interaction. 



The present study therefore provides a detailed base- 
line calculation for the search of critical phenomena 
in event-by-event multiplicity fluctuations. If non- 
monotonous deviations from these predictions, as e.g. 
expected by droplet formation, are observed, these en- 
hanced fluctuations might indicate the onset of decon- 
finement and/or the critical point. 

The NA49 collaboration is currently studying the 
energy dependence of multiplicity fluctuations from 
20 A GeV to 158 A GeV . Further detailed exploration 
is planned for the NA61 (SHINE) experiment 51] at the 
CERN SPS. In addition the critRHIC experiment and 
the CBM experiment [zl at the GSI FAIR facility as 
well as and the MPD experiment at JINR in Dubna will 
be able to explore this energy region in the near future. 
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